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1 Introduction 

Let {X n } n >i be a stationary sequence with marginal distribution Fx n = F. 
Consider M l ;j = \/ J s=i+1 X s , where x\/ y denotes max(x,y), with Mqj = Mj 
and Mij = — oo for i > j. The sequence {X„}„>i has extremal index 9 G [0,1] 
if, for each r > 0, there is a sequence of normalized levels {u n = ui T ^}„>i, i.e., 

n(l - F(u n )) -A r, 


as n —> oo, such that 

P(M n < u n ) —> e 6t , (1) 

(Leadbetter et al. [18] 1983). When 9 = 1 the exceedances of high thresh¬ 
olds u n , by the variables in {X n } n >i, tend to occur isolated as in an inde¬ 
pendent variables context. However, if 9 < 1 we have groups of exceedances 
in the limit. Clusters of extreme values are linked with incidences and dura¬ 
tions of catastrophic phenomena, an important issue in areas like environment, 
finance, insurance, engineering, among others. The extremal index is a key 
parameter in this context and its estimation has been greatly addressed in lit¬ 
erature. The most popular procedures are the blocks and the runs estimators 
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(e.g., Nandagopalan [JTj 1990; Hsing |12| 1993; Weissman and Novak |28| 1998; 
Robert et al. 22 2009), and more recently, the interexceedance times method 
(Ferro and Segers [TO] 2003; Siiveges [27|). The first group estimators requires a 
clustering identification parameter which is a largely arbitrarily task to comply 
and has some impact in inference. The second group avoids this parameter but 
may be less precise in specific contexts. 

In this paper we propose a new estimation procedure that works under the 
local dependence condition D^ k \u n ) of Chernick et al. ([0], 1991). This condi¬ 
tion requires the dependence condition D(u n ) of Leadbetter ([IT], 1974), which 
states that a n ,i n —> 0, as n —> oo, for some sequence l n = o(n), where 

&n,L — SUp{ | P ( ^ — ^n) P (i ,ii+p Un) P ~ ^n) 

1 < h < k + P + l < ji < ji + q < n}. 

We say that condition D( fc )(u„) holds for {I„}„>i, if for some {fc ra }n>i such 
that, 


k n ^ oo, k n otn,i n ^ 0, k n l n /n ^ 0, 

as n —> oo, we have 


( 2 ) 


nP (Xi > u n , Mi jfc <u n < M k rn ) —0, 

n—>■ oo 

with { r n = [n/k n ]} n >i ([x] denotes the integer part of x). Condition D( fc )(it„) 
is implied by 


r n 

n P (Xi > u n , M lk <u n < Xj) —0. 

z — J n—> oo 

j=k +1 


This corresponds to condition D ( u n ) of Leadbetter et al. f[18|. 1983) whenever 
k = 1 which locally restricts the occurrence of clusters of exceedances and 
thus leads to 0 = 1. If k = 2 we have condition D ( u n ) of Leadbetter and 
Nandagopalan ([IS], 1989). This condition locally restricts the occurrence of 
two or more upcrossings, but still allows clustering of exceedances. 

In Chernick et al. ([3], 1991) it is proved that, under D( fe )(u n ), the extremal 
index exists and is given by 

0 X = lim P{M l k < u n \Xx > u n ). (3) 

n—>oo 

The runs estimator can be derived from this relation by taking the runs param¬ 
eter r as k. In particular, under condition D( 2 )(wn)> Nadagopalan (22, 1990) 
found 

0 X = lim P(X 2 < u n |Xi > u n ) 


P(Xi < Un < x 2 ) 

= lim --- j -, 

n->-oo P(X 2 > U n ) 

which motivates his estimator based on the ratio between the number of upcross¬ 
ings (equal to the number of downcrossings) and the number of exceedances. 
Although the D^^rin) condition implies D( fc )( M n) for k > 2 and we have several 
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representations for 9x as in ([3]), under D^ 2 ^(u n ) we have only to be concerned 
with the count of upcrossings and exceedances, rather than the length r for runs 
of non-exceedances or intervals between exceedances. It is this easy approach in 
the Nandagopalan’s estimator that we want to take advantage in this paper, by 
estimating 9x through the extremal index of an auxiliary sequence satisfying 
D (2) (Urc)- 

The results that motivate our new estimation approach are given in Section 
[2j In this section we relate the extremal index 9x of the process satisfying 
D ( k \u n ) with the one of a process of cycles satisfying D^ 2 ) (u n ), deriving new 
representations for 6 that motivate the estimators. In this way, we promote the 
application of the estimation procedures that work under D ^ 2 \u n ). 

Knowledge about D^ fc ^(u„) has not only impact on the computation of the 
extremal index of a process but also informs about the cluster structure of ex¬ 
treme values. In moving maximum processes we directly obtain the extremal 
index by calculating the limit in Jl]). It may be the reason why there is no study 
in literature, as far as we know, concerning local dependence within these pro¬ 
cesses. In Section [3] we derive a necessary and sufficient condition for (u n ) 
to hold within moving maxima processes. 

Section 0] is devoted to inference. We state a diagnostic tool to analyze 
( u n ) since it is the context of our framework. Therefore, we are also moving 
forward in diminishing the arbitrarity in the declustering scheme of the runs 
estimator. We analyze the performance of the new estimators trough simulation 
and illustrate with an application to a financial time series. We conclude in 
Section [5] 


2 Extremal index of grouped variables 

Let {/o = 0,/ n ,n > 1} be an increasing sequence of integer random variables 
(r.v.’s) such that {S'™ = I n — I n -i}n>i is an i.i.d. sequence satisfying E(S n ) = p. 
From such a renewal process and a stationary sequence {X n } n >\, let define 

Z n = M In _ lJn , n> 1. (4) 

Driven by the strategy used by Rootzen (1988, j23]) in the study of the extremal 
behavior of the regenerative processes, we will compare M n with the maximum 
of the first [n/p] variables in the sequence of cycles {Z n } n > i. 

Proposition 2.1. Let {X„}„>i be a stationary sequence and {Z n } n > \ de¬ 
fined by for some renewal process {S , n }n>i such that E(S n ) = p. If 

l n \li >i P {Mi i _ 1 ,i i > Un) r is bounded, then 

{ — J n>l 


P(M n <u n ) - P 




—► 0, n —> oo. 


Proof. Let L n = sup{fc : Ik < n} and U n = infjfc : Ik > n}. By the law of large 
numbers, Ve > 0, we have 


( 

Ln_ 1 

> e ] ->• 0 and P ( 

V 

n p 

J n—t oo V 



0 . 


P 


n—> oo 
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Furthermore, 

p(M n < Un ) <p (n,i"i{Af w . < “«}) 

= p (r\f=i{ M h-i,ii < U n },n(l/p- e) <L n < n(l/p + e)) + o(l) 

<p(nWM{M W( <« B })+o(i) 

Similarly, we derive 

p(M n < Un ) > p (ni=i 1/1H ' e)1 {Mj 1 _ 1 ,j i < «„})+o(i). 

Now, just observe that 

o < p (ctJihMu^u < Un }) - p (n!: ( i 1/p+e)1 {M /l _ 1 ,/ i < «„}) 

< ne \li>\{P {Mii_iji >u n )} < ek, 
as well as 

o < p (ni= ( i 1/p " e)] {Mz i _ 1 ,i i < Un }) - p (ni: / 1 Pl {^ i - 1 ,/ i < u n }) < ek , 

for some constant k. □ 


If we assume that {Z n } n > i is stationary satisfying D(u ra ) and a local de¬ 
pendence condition D( fc )(M„) then, by applying the previous proposition, we 
can compute the extremal index of {X n } n >i from the knowledge of the joint 
distribution of a finite number of consecutive terms of {Z n } n > i. 

In what concerns the local behavior of the large values of {Z n } n > i, we are 
going to consider to ways: in Proposition 12.21 we derive the extremal index by 
assuming the local independence condition DW(u n ) and in Proposition 12.31 by 
assuming the local dependence condition D( 2 )(it„). 

Proposition 2.2. Under the conditions of Proposition\Kl} if {Z n }n>\ is sta¬ 
tionary and satisfies D(u n ) and D W (u n ) conditions for u n such that un = 
for {X n } n >i and u n = uff ^ for {Z n } n >then {X n } n >i has extremal index 

e x = iim > = IL. 

n^oo pP(X i > U n ) PT 

Proof. We have that {Z n } n > i has extremal index Oz = 1 and thus, by applying 
Proposition 12.11 


lim P(M n < u n ) 

n—>• oo 


lim P 

n—>• oo 



lim (P(M h < u n )) 

n—> oo 


W/p] _ e ~ T * /P' 


Therefore, lim P(M n < u n ) 
Ox 


= e SxT , as n —>• oo, with 
t* P{Mi i > u n ) 

= — = lim —--—. 

rp n->oo P[X 1 > U n )p 


□ 
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This is what happens in regenerative processes with independent cycles (see 
expression (4.2) in Rootzen, [231 1988, obtained directly). 

Proposition 2.3. Under the conditions of Provosition WU\ if {Z n } n > i is sta¬ 
tionary and satisfies D(u n ) and D^ (u n ) conditions for u n such that un = uff ^ 
for {X n } n > i, u n = u[f } for {Z n jn>l a,nd nP (M h <u n < M/ 1i/2 ) ->• v*, then 
{X n } n >i has extremal index 

a 0 Z t* P(M h < u n < M Ilj2 ) v* 

U x =-= inn -----= —. 

pT TI-HX) pP(X 1 > u n ) pr 


Proof. We have that {Z n } n > i has extremal index 

P{M Il < u n < M Ii j 2 ) v* 

uz = hm - „., r -r-= — 

n-too P(M/ 1 > U n ) T* 

and thus, by applying Proposition l2.il 


[n/p] 


lim P(M n < u n ) = lim P I \J Zi < 

n. — n .—^no l ' 


Un I = e-Oz^/p = e -re x 


with 


Q OzT P(Mi 1 < Un < 

OX = - = llm -7TTT7-.-■ 

pr n-xxs pP(Xi > U n ) 


(5) 


( 6 ) 

□ 


This is what happens in regenerative processes with 1—dependent cycles, 
(see comment after expression (4.2) in Rootzen, (233 1988). 


Since 


TlP I Z\ Un, } \J Zi ~^> Un 


i -2 


= nP 


^^1 ^ U n , Z 2 U n 


<\J zA +nP (Z 1 > U n , Z 2 > U n ) 

i—3 ) 


and 

nP ( Z 1 > u n , Z 2 > u n ) = nP (Z x > u n ) - nP {Z 1 > u n > Z 2 ), 

we can remark that, for {Z n } n >i satisfying D < ' 2 )(u rl ), it holds that { Z n } n >i 
satisfies D^^zi^) if and only if r* = v*, that is, the limiting mean number of 
exceedances is asymptotically equal to the limiting mean number of upcross- 
ings (or downcrossings). Also, for any k > 2, provided that fX n } n >i satisfies 
D (k) {u n ), it holds (u„) if and only if 

lim n ( P (Ad > u n , Mi fc _i <u n ) - P (Ad > u n , Mi jfc < u n )) = 0. 

n—> 00 

This remark will help in the choice of a value for k, in Section [H dedicated to 
the estimation of Ox- 
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The presented results also point out a way to obtain the limiting law of the 
maximum term of the first ^" =1 Si r - v -’ s of the sequence {X„} n >i. In fact, for 
{X n } n >i and {S'n}ra>i as above, it holds that 


/£r=i* \ 

P 1 Si < Wn) = I V Xj <U n \= P(Zi < u n ,..., Z n < u n ) 

—)• e~ 9zT = e~ pT8x = exp {— E(S\) lim nP(X i > u n ,X 2 < u n ,... ,Xk < u n ) 

—^oo l n—>oo 


if {X„}„>i satisfies D ^ k \u n ). 


We could state a general result analogous to the above Propositions 12.21 and 
12.31 by considering {Z n } n > i satisfying with k > 2. However, our final 

goal is to relate the extremal index of a sequence {X n } n >i satisfying D ( k \u n ) 
with k > 2, with the extremal index of an auxiliary sequence {Z n } n > i satisfying 
D^(u n ) with k < 2. This will enable to take profit of the estimation of the 
extremal index under or D^ 2 ^(u n ), after a suitable transformation of 

the data. The identification of clusters reduces then to the identification of 
blocks of consecutive exceedances. Next results discuss relations on long-range 
and local dependence conditions for {X n } n >i and {Z n } n > i, which can be easily 
obtained in the particular case of a deterministic n > 0, considered later for 
the main proposal of this work. 


Proposition 2.4. Let {X n } n >i be such that { Xi , i £ B} is independent of 
{Si, i £ A} whenever An B ^ 0, for some renewal process S = {/o = 0, S n = 
In — fy_i} n >i with t < S n < s, n > 1. If, conditionally on S, {X n } n >i satisfies 
condition D(u n ) with spacer sequence l n , then {Z n j n >i, defined by W’ satisfies 
D(u n ) with 1* = [2l n /t\. 


Proof. Let / = {ii,...,i p } and J = {ji, ■ ■ ■ ,j q } with 1 < fy < ... < i p < 
i p + Z* < ji < ... < j q < n and l* = [2l n /f\. Then 

| P (Vie/ u m Viej S Un) — P (Vie/ %i — u n ) P (VigJ < U n ) | 

= P (j 3 (Vie/*(Sj) Xi < u n , Viej*(Sj) Xi < u n ) |<S^ 

~P (P (VieiASA Xi < u n ) |s) E (p (V ieJ . (Sj) X t < u n ) |s) , 

where Sa denotes the vector of r.v.’s Si, i £ A, and I*(Si) and J*(Sj) are 
separated by at least l n , since we have Jj < i p s, Ij 1 -i > (ji — 1 )t and thus 
(ji — l)t + 1 — i p s > l n , for large enough n. Therefore, and meeting the given 
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assumptions, the previous expression is upper bounded by 
E (p (y^ (Sl) Xi<u n ,\j i£j* (Sj) X i ^ U n ) \Siujj 


-E 
+ E 
-E 

&n,L 


(j 3 (Vie/* 


(Sj) Xi — 


(Sj) - u? 


(Vig/*(Sj) <U n ) P 


(^(V i£ z. 


(S/) X i — U 1 


)p(yiej* 

(v, eJ . 

) \Si) e(p (ViGJ* ( s 


(Sj) X i — U * 


) I-S'/uj) 


j j Xi ' ■ U r 


□ 

Proposition 2.5. Let {X„} n >i fee a stationary sequence and {Z n } n >\ defined 
by I^J), for some renewal process S = {Iq = 0, S n = I n — / n -i}n>i with t < 
S n < s, n > 1. 

(a) If {Z n } n > i satisfies (u n ) with r n = [n/fc„], t/ien {}n>i satisfies 

D (2s-t+l) 

(u n ) with the same r n . 

(b) If 2t > s and {X n } n >i satisfies (u n ) for some k < 2t — s + 1, with 
r n = [n/k n \, then {Z n } n > 1 satisfies D (2) (u n ) with r* = [r n /s\. 

Proof. To obtain (a) we take into account the following inequalities: 

TlP {X i > U n , Mi^s-t+l < U n < M2s-t+X,r n ) 

< nP ( X\ > U n , Mi^s-t+l f U n < M2s-t+l,tr n -t+l)) 

= TlP (Xt U n j Mt,2s if n n <1 Xl2s,tr n ') 

< nP (Mi 1 > u n ,M Ix j 2 <u n < V [=3 Mii-iji) 

= o(l), n —> oo . 

For (b) we have: 

nP ^Mi x > u ni M Ix j 2 <u n < V [=3 Mfi-iJi'j 

< nP (^M s > U n , M s , 2 t < U n < Vi= 2 t +1 X j) 

< n X^j=i P (Xj > n n , Mj s < u n ,M Sy 2t <u n < M 2 t,r n ) 

and each of the s terms in the sum above tends to zero by the D^ 2t_s+1 )(u rl ) 
condition for {X n } n >i. □ 

We state now a result on the “clustered" process {Z n } n > i that resumes the 
path to obtain its extremal index 9z and to recover 9x for the “declustered" 
process {X n } n >i, enhancing that to count the mean number of upcrossings (or 
downcrossings) for {Z\,..., Z n } is asymptotically equivalent to count the mean 
number of runs {Xi > u n , X i+1 X l+fe _i < u n }, i < np. 






Corollary 2.6. Let {W„}„>i and S be in the conditions of Propositionand 
that E(S n ) = p. Suppose that u n = Un' 1 for {X n } n >i and {X n } n >i satisfies 
D^ (u n ) for some k < 2t — s + 1. Then 

(a) {Z n } 7 !>i defined by (0) satisfies D(u n ) and (u n ) conditions. 

(b) Ifu n = Un ■* for {Z n } n >i and there exists lim nP (Mi 1 <u n < Mi ± j 2 ) = 

n—too 


v* then it holds that 


and 


_»* _°ZT* 

VZ — —, Vx — - 

T* pr 


lim nP (Mj 1 < u n < M h j 2 ) = p lim nP (Xx > u ni X 2 < u n ,..., X k < u n ). 


We now focus on the particular case of I n = n(k — 1), n > 0, for some k > 2. 
Therefore, we have s = t = k — 1. If is this the case then, from the previous 
result, condition D^ 2 ^(u n ) for {Z n } n > i implies condition D ^ k \u n ) for {X n } n >i 
with the same length r n to model “local behavior". Otherwise, the validity of 
condition D ^ k \u n ) for {X n } n >i leads to condition D^ 2 ^(u n ) for {Z n } n > i with 
r* = [r n /[k — 1)]. We can then state the following corollary, which can also be 
proved directly. 


Corollary 2.7. Let {X n } n >i be a stationary sequence and {Z n } n >i defined by 
Op with I n = n(k — 1), n > 0, for some k > 2. Then {Z n } n >i satisfies condition 
L F 2 ) (u n ) if and only if {X n } n >i satisfies condition D^ (u n ). 

Thus, under conditions of Proposition ^. 5l and according to (JG]) , the extremal 
index can be written as 


Ox 


OZT* 

(k — l)r 


= lim 

n—> oo 


PjMl! < u n < M hj2 ) 
(k - 1 )P(X 1 > U n ) 


(7) 


lim 

n—too 


P(Mk -1 <u n < M k _i t 2 (k~ i)) 
(k - 1 )P{X 1 > Un) 


By using the stationarity, this limit can be rewritten as the one obtained in 
Chernick et al. ( 0 , 1991) under condition and given in (J3J). However, 

representation 0 allows to estimate 8x through the estimation of 9z for the 
cycles Zi = i > 1, for which D^ 2 )(u n ) holds, as we will present in 

Section |U 


Here we illustrate the above results with finite moving maxima processes 
and in the next section we devote special attention to the local dependence in 
this kind of processes. 

Example 2.1. Consider the moving maximum process, X n = \f )=() 
ao = 2/6, ai = 1/6 and a 2 = 3/6, n > 1, with sequence {Y n } n >-1 independent 
and having standard Frechet marginal distribution, Fy = exp(~l/a;), x > 0. 
This stationary sequence satisfies D^ 3 )(u n ) for levels u n = n/r, r > 0, as will be 
seen in the next section and 9x = ((2/6) V (1/6) V (3/6)) = 1/2 (see Weissman 
and Cohen 05], 1995). For Z n = X 2n ~i V X 2n , we have nP{Z\ > u n ) —> 
(5/3)r = r*, as n —>• oo, since 

P(Z 1 < u n ) = F Y (u n )F Y (3u n /2). 
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Observe also that 

P(Z 2 < u n ,Z 1 < u n ) = F Y (u n /2)F Y (3u n /2) 
and, provided that D^(ii n ) holds for {Z n } n > i, we have 

0 Z = lim P(Z 2 < u n \ Z\ > u n ) = 3/5. 

n—»oo 

By applying 0 , we obtain 9x = 1/2. □ 


We have seen that, for every stationary sequence {X n } n > i satisfying (u„), 

k > 2, we can build a stationary sequence {Z n } n >i satisfying D^ 2 )(u n ) by tak¬ 
ing the maxima of k — 1 consecutive variables of sequence {X n } n >i. For big 
values of fc, such aggregation can result in reduced accuracy in the estimation 
of 9x via the sample based on {Z n } n > i, as will be pointed in Section [4] The 
Proposition 12.51 also states that it can be considered mixtures of big cycles of 
several lengths in order to build the sequence {Z n } n >i satisfying D^ 2 )(«,„), as 
it is illustrated in the next example. 

Example 2.2. Let S = {Iq = 0,5 n = I n — I n -i}n>i be a sequence of in¬ 
dependent r.v.’s uniformly distributed on {k, k + 2}, with a fixed k > 6, and 
independent of {X„}„>i. For {X n } n >\ take a stationary sequence satisfying 
D( 5 )(w! t ' ) ) and D(ui T ' ) ), for instance, a MM process with signatures aij as given 
in the next section and u = n/r. Let Z n be as in (U, which is stationary and 
also satisfies D(u n ). Since 2k — (k + 2) +1 > 5, then {Z n } n >i satisfies condition 
D( 2 )(u n ). It holds that 


T 


* 


lim — 

n—too 2 


nP 

s€{k,k-\- 2 } 



v* = lim - V 

n—too 4 J 

si,S2€{k,k+2} 


nP 


\/*i< 



Oz = yv and 9 X 


IS* 

(fe+l)r ' 


□ 


3 Condition (u n ) for max-stable processes 

A moving maxima process (MM) is defined as 

X n — \J V a kj Y i,n-j, n> 1 (8) 

/>! — oo<j <oo 


where {Yjj, l > 1, —oo < j < oo} is an i.i.d. sequence of r.v.’s, usually unit 
Frechet and {cxij, l > 1, — oo < j < oo} are non negative constants (usually 
denoted signatures) such that J2i>iJ2-oo<j<oo a i,j = 1 (Deheuvels [Bj 1983, 
Davis and Resnick [5] 1989, Smith and Weissman |20l 1996, Hall et al. [II] 
2002, Meinguet [2D] 2012). 
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Under the condition ]C-oo<j<oo Sz>i the MM process is strong¬ 

mixing (Meinguet, [2D] 2012) and therefore it satisfies D(u n ). 

An interesting feature of these processes is that the transformation of {Yi j, l > 
1, —oo < j < oo} induces a dependence structure with extremes in temporal 
clusters. Any stationary process with finite distributions of multivariate extreme 
value type can be approximated by an MM process with marginals of extreme 
value type (Hall et al. [TT] 2002). Examples of finite MM processes (i.e., with l 
and j finite) are not difficult to deal with and are often used to illustrate long 
range and local dependence conditions within extreme values. The extremal 
index is directly obtained through lim^^oo P(M n < n/r ), r > 0, even for infi¬ 
nite MM and thus avoid the validity of some condition. In Meinguet ([20] 
2012) it was presented a nice finite-cluster condition which prevents a sequence 
of extremes occurring in MM from being infinite over time. However, it doesn’t 
enable a representation for 6 from finite marginal distributions of the process. 
The local dependence conditions brings us enlightenment about the clustering 
structure of extreme values. Any finite MM is m-dependent for some positive 
integer to and thus holds at least for some k > to. From simple examples, 
we know that small changes in the values of coefficients aij may lead to large 
diferences within the clusters structure. Hence, it is raised the question of which 
conditions must satisfy cqj so that some holds for an MM process. The 
next result presents a necessary and sufficient condition. 

Proposition 3.1. Let {X n } n >i be an MM process as defined in (0), where 
l > 1, —oo < j < oo} is an i.i.d. sequence of unit Frechet r.v. ’s. Then 

(a) {X n } n >i satisfies condition D^(u n ), k > 2, if and only if, for all l > 1 
and — oo < j < oo, 

a l,j A I \J i J < \J oti^j+ s - 1, (9) 

ys>/c+l J s —2 

where x Ay denotes min(a;, y ). 

(b) {X n } n >i satisfies condition (u n ) if and only if, for all l > 1 there is 
only one — oo < j < oo such that cqj > 0. 

Proof. (a) Observe that 


nP I A i > u n , X s f u n , X s > Uyj 


s=2 


s=/c+l 


= nP Al > u n , V X s > u n ) - nP ix 1 > u n , \JX s >u n , V X s > U n J . 


s—k -\-1 


s —/c—|— 1 


Since X s = Vz Vj a l,jYi,i-(j- s+ i) and {Yij, l > 1, -oo < j < oo} is a 
sequence of independent r.v.’s, we have 


lim nP Ai > u n , \J X s > u r 

n—± oo \ v 

\ s—k-\-l 


( 10 ) 


= lim n > > P 

n —Voo • ^ ^ 


i j 


\f &l,j+s—l^ 1,1—j 

s=k-\-1 




11 


Thus, for levels u n = n/r, t > 0, 

lim nP I Xi > n/r , \f X s > n/i 

n.—¥ oo \ V 


s=k+l 


l 3 


= lim n Y Y P (Yi iV 

n .—loo « V 5 /'v < • 


V 7 


= i n EE p y, 'H> 




Vs=fc+l QI| J+ 8_1 

n/r 


a i,3 A V s =fc+1 a i,i+s-i / 


Using the same reasoning on the second term in (1101) and by the theorem 
of dominated convergence, 

/ k r„ \ 

lim nP I Xi > u n , V 


s—2 


s=fc+l 


= 

n/r 


n/r 


i j 
— Fy 


a l,j A Vs =2 a iJ+s—1 A Vs=fc+l a i>J+s— 1 . 


a hl A Vs=fc+1 a 0i+s-i 


= y E T aij A V a i j+s-1 - a U A V “hPM-i A V a bPH>-i 


l j 


s>k +1 


s>fc+l 


since 1 0, as l —> oo and j —1 oo. Now just observe that the limit will 
be null whenever relation © holds, for alll > 1 and — oo < j < oo. 

(b) The conclusion in (b) follows from 


lim nP I X i > Un, \J X s > u n 

n—>• oo \ v 

\ 8 =2 / 


EE r ( aij A V COj'+s-l 

i j 


s >2 


The latter sum is null if and only if, for each l , there is j* such that 
a;j. > 0 and aij = 0 for j ^ j*. 

□ 

Example 3.1. Consider the moving maximum processes, X n = \l^_ 0 ajY n _j, 
ao = 2/6, ai = 1/6 and «2 = 3/6, n > 1, given in Example 12.21 and 
W„ = \j3_ 0 ajY n _j, ao = 1/6, a\ = 3/6 and a .2 = 2/6, n > 1, with sequence 
{Y n } n >- 1 independent and having standard Frechet marginal distribution. We 
will see that {X n } n >\ satisfies D^ 3 ^(it„) (and not D^ 2 ^(m„)) and that {W n } n >i 
satisfies D^ 2 ) (u n ), for levels u n = n/r, t > 0, by applying relation ©. The 
calculations are summarized in Table[T| Observe that D( 2 )(u„) doesn’t hold for 

{X n } n >i since, if j = 0 then ag A (V s >3 a s-i) = 2/6 A 3/6 > 1/6 = a.\. □ 
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Table 1: Verification of conditions D^ 3 ^(u„) and D( 2 )(u n ) for, respectively, the 
MM processes {X n } n >i and {W n } n >i of Example 13. 11 according to relation in 

ED- 



condition D( 3 )(it n ) 

for 


condition D ^(u n ) 

for {W n } n >i 

j 

a j A (\/s>4 

otj+ 1 V aj-t-2 

3 

a j ^ ( Vs>3 a j+s-l ) 

a i+ 1 

< -3 

0 A 3/6 

0 

1 

VI 

0 A 3/6 

0 

-2,-1 

0 A 3/6 

2/6 

-1 

0 A 3/6 

1/6 

0 

2/6 A0 

3/6 

0 

1/6 A 2/6 

3/6 

1 

1/6 A 0 

3/6 

1 

3/6 A0 

2/6 

2 

3/6 A0 

0 

2 

2/6 A0 

0 

> 3 

0 A 0 

0 

> 3 

0 A0 

0 


Inference within MM processes has been addressed in literature (see Zhang 
and Smith [3D], 2010). Therefore, as an alternative to the empirical method of 
Siiveges (EH 2007), we can check the validity of D( fc ) (u n ) within these processes 
by estimating coefficients aij and applying ED- 

The MM processes are stationary max-stable processes for which, under 
D( fe )(un), we can derive the extremal index from a tail dependence coefficient. 
Suppose that the stationary process { X n } n >i has unit Frechet marginals F(x) = 
exp(— 1/x), x > 0. If D( fc )(u^) holds for {X n } n >i, then 

Ox = lim P (Mi k < n/r\Xi > n/r) 

= 1- lim P (Mi k > n/r\Xi > n/r) (11) 

n—*oo 

= 1-A^(1,1), 

provided the limit exists, where I\ = {2,..., fc}, I 2 = {1} and A'| l|/2) (1,1) is the 
upper tail dependence coefficient considered in Ferreira and Ferreira ([S], 2012b). 
In the case of max-stable processes or, more generally, processes satisfying the 
max-domain of attraction condition, the limit in m is always defined. By 
applying the propositions 2.1 and 3.1 in Ferreira and Ferreira ([D], 2012b), we 
conclude that 


E[e~ M ^ E(e~ M ^ 

X 1 -E 1 - E 

and, in particular for k = 2, it holds that 

° X ~ 1 - E (e-( x ' vx z')- 1 ) ~ 2 


( 12 ) 


(13) 


This representation for Ox motivates its estimation from moment estimators 
iov E = E (yt.FiXi)), as considered in Ferreira and Ferreira ( 0 , 

2012b). 
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We apply now the results of the previous section in order to compute 6x 
from d z of the process {Z n = V"i(n-i)(fe-i)+i X i}n>i under D (2 >(w„). The 
estimation of 0 Z is considerably more simpler as suggested by m 

Proposition 3.2. Let {X n } n >i be a stationary max-stable process with unit 
Frechet marginals F and u^ = n/r, t > 0. Then 

(a) {Z n } n >i is stationary and max-stable with marginal distribution F z {x) = 
F ek ~ 1 (x), where e k -i = - logi i (x 1 ,... > x fc _ 1 )(l, • • ■, 1) £ [1, A; — 1] is the 
(k — 1 )-th extremal coefficient of {X n } n > i. 

(b) If {X n } n >i satisfies D(u n ) and D^(u n ), k > 2, then {Z n } n >\ satisfies 
D(u n ) and D ^ (u n ), 


1 -E(F z (Z 1 )\/F z (Z 2 )) 

and 

-logFz(l) 
dx ~ 0Z k-1 ■ 


(14) 


(15) 


Proof. We only justify (b), leaving (a) to the reader. 

We first consider the sequence of cycles {Z* = Vri(n-i)(fc-i)+i x i/ e k-i}n>i 
which satisfies the same local and long-range dependence conditions as {Z n } n > i. 
For this stationary and max-stable sequence with unit Frechet marginals, by 
applying m we obtain 

° Z ‘ = l ~ E ~ 2 ' 

Then 

dz = ° Z * = 1 — E (e - ( ZlVZ2 ) -l£fe -i) “ 2 = 1-E{F z (Z 1 )VF z {Z 2 )) ~ 2 ' 

To obtain the relation (fl5l) we apply Proposition 12.31 with 
r* = lim nP{Z > n/r) = —\ogF z (l)T. 


□ 

This result suggests the estimation of Ox via the estimation of — logPz(l) 
and E ( F z {Z\ ) V F z (Z 2 )). 

4 Estimation 

Our new estimation proposal consists in first, to state the sequence of cycles, 
Z n = VeL k (n-i)(k-i)+i X s ’ n — anc ^ tl ien estimate 0 based on {Z n } n > 
Observe that, from Proposition 12.31 we can define the estimator 

n = U nM 

x NX(u n y 


(16) 
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as well as, the estimator 


e z N%(u n ) 

N*{u n ) 


(17) 


where U^{u n ) and N^{u n ) are, respectively, the number of upcrossings of u n 
and the number of exceedances of u n within {Z i,..., Z\ n i(k-x)\\ and N^(u n ) 
is the number of exceedances of u n within { X \,... ,X n }. Since D( 2 )(tt n ) holds 
for {Z n } n > i, estimators under this condition can be used to calculate 9 Z , e.g., 
the maximum likelihood in Siiveges ( EZ] , 2007) and the upcrossings estimator 
in Nandagopalan ([23], 1990). 

Now observe that, based on ED. we can write 0z as 


0 Z = 1 — lim 

n—> oo 


> U n \M Il > U n ) 
P(M/ 1 > U n ) 


= 1 — A z , 


where X z is the so called “tail dependence coefficient" (see Joe [H] 1997 p. 33, 
Cotes et al. |4j 1999, Schmidt and Stadmiiller [24] 2006 and references therein; 
see also Ferreira and Ferreira [5] 2012a Proposition 4). Hence, we can derive 


e x 


(1 - A z )t* 
(k — 1 )r 


and thus also state the estimator 

2 (1 ~ Xz)NZ(u n ) 

X NX(u n ) 


(18) 


(19) 


We can estimate the tail dependence coefficient by applying a non-parametric 
procedure, e.g., the one in Schmidt and Stadmiiller ([23], 2006). For the par¬ 
ticular case of max-stable processes, by representation ED, we can also apply 
the estimators 9 Z in Ferreira and Ferreira ([jj], 2012b) for 9 Z and again Ox 
as in ED- A similar procedure based on (fT5l) leads to a second estimator for 
max-stable processes, namely 


2 2 -logPz(l) 

9x ~ 9z k- 1 ’ 

where F Z ( 1) is the empirical distribution function. These estimators will be 
denoted, respectively, 9 ss , 6 FF and 9 FF . 

In the next section we analyze our new proposal through simulation. For 9 Z 
in expression ED, we consider the upcrossings estimator of Nandagopalan (eh, 
1990), the estimator of Ferro and Segers ([ID], 2003) also known as intervals es¬ 
timator and the maximum likelihood estimator of Siiveges (EZJ, 2007), and 
denote our extremal index estimators, respectively, 9 17 , 9 1 and 9 ML . We also 
compare with the intervals and runs estimators applied directly on {Al„}„>i. 
For these estimators we use notation 9 1 and 9 R , respectively. 

In order to analyze D( fc ((«.„) and construct the cycles {Z n } n > i, we can ex¬ 
tend the methodology in Siiveges ([27], 2007) considered to check D( 2 )(u n ). 
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More precisely, we compute the proportion of anti-D( fc ) (u n ) events by 

E n— r„+l 

7 — 1 f A.J , Aj-1-1 . .,Xj_ j_fc — x <.U n , Mj+k _1 r —1^77^} 

Pk[u n ,r n ) = -=— -, 

Z_7j=l 1{A- 3 >u„} 

for normalized levels u n approximated by the empirical quantiles 1 — r/n, 
for some fixed positive r and some sequence {r n } n > i satisfying the condi¬ 
tions of Proposition (12.21) . We take the proportions Pfc(it m ,r m ) for sequences 
{Xi,, X m }, with increasing length m < n. 

On what concerns the choice of {k n } n > i, we can choose, for instance, the 
family of sequences of integers kn' 1 = [(logn) s ], s > 0. Thus, for each r and s, 
we can plot the points (m,pfe(u m , rm )), which must converge to zero, for some 
s, as m increases if (u n ) holds with k„^ ■ This is a slightly different approach 
of the one in Suveges ([27], 2007), but closer to the definition of («.„), since 
this condition states a limiting behavior as n -> oo, and u n « ,F _1 (1 — r/n), 
r n = [n/k n \ and p k {u n ,r n ) are functions of n. To avoid three-dimensional plots 
that arise from the joint variation of r, s and m, we can separately analyze the 
evolution of the proportions for different choices of r n . For the particular case 
of (un) condition, we have the proportions 


Pi(u n ,r n ) 


E n—r n +1 -f| 

j— 1 \>Un — 1 } 


E Tl 
7 — 


3 =1 


1 ( *j>u n } 


Once accepted the condition ( u n ) for some fc, that means we consider that 

the process satisfies D^ s )(« n ) for all s > k and does not satisfy D( s )(it„) for 
s < k. The decision to exclude values less than k may be based on the analysis 
of (rn,pk-i{u m ,r m )) or, from the remark after Proposition 12.31 by comparing 
d k -i{u n ,r n ) with d k {u n ,r n ), where 


n-r n -\-1 

d k (u n ,r n ) = l{Aj>ti n ,Mj, J+fc _i<u n }' 

j=i 

A good choice of k is enhanced by away trajectories for (m, dk-i(u m , r m )) and 
(m,d fc (zi m ,r m )) and close trajectories for (m, dk(u m ,r m )) and ( m,dk+i{u m ,r m )). 

An illustration is given in Figures [1] and [2] where it was considered a simu¬ 
lated sample of size 10000 for each of the models given below. In Figure [l] it is 
plotted the proportions of anti-D^ 3 ) ( u n ) (first five panels) and anti-D^ 4 ) ( u n ) (last 
panel), with k n = [(logn) 3 ] and values r = 50,100, for some models. More pre¬ 
cisely, the three first line panels correspond to the proportions of anti-D® ( u n ) of 
the following models: a first order autoregressive process with Cauchy marginals 
and autoregressive parameter p = —0.6 of Chernick ([ 23 , 1978), a negatively 
correlated uniform AR(1) process of Chernick et al. ([3], 1991) with r = 2, 
respectively denoted ARCauchy and ARUnif, and an MM process with coeffi¬ 
cients «o = 2/6, a\ = 1/6, «2 = 3/6 as given in Example 12.21 The first second 
line panel correspond to the proportions of anti-D^^itn) of a first order MAR 
process with standard Frechet marginals and autoregressive parameter <f> = 0.5 
of Davis and Resnick ([5], 1989). The last two panels in the second line corre¬ 
spond to, respectively, the proportions of anti-D^ 3 ) (it„) and anti-D^ 4 )(it„) of a 
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Markov chain with standard Gumbel marginals and logistic joint distribution 
with dependence parameter a = 0.5 The MAR process satisfy D^ 2 )(u„) and so 
D ( 3 >(u n ) holds, thus leading to small proportions of anti-D( 3 )(u n ). The same 
scenario is noticed in the first three models, all satisfying condition D( 3 )(u„). 
There are slightly upper curves within the Markov chain but still comprising 
small values. A little decrease occurs in the proportions of anti-D^ 4 )(M n ) for the 
Markov process. 

In Figure [2] we find the proportions of anti-D^ 3 ) (u n ) to anti-D^ 5 )(u n ) of a 
GARCH(1,1) process with Gaussian innovations, autoregressive parameter A = 
0.25 and variance parameter /? = 0.7 (Laurini and Tawn, [1(3] 2012). More 
precisely, in the first two panels are plotted the proportions of anti-D* 3 ^ (u n ) 
by choosing k n = [(logn) 3 ] and k n = [(logn) 3 ' 3 ], respectively, and the last two 
plots correspond to the proportions of anti-D*- 4 ) (u n ) and anti-D'- 5 ) (u n ). with 
k n = [(logn) 3,3 ]. We can see that the choice k n = [(logn) 3 ' 3 ] may be better 
within this case. The plots also suggest that condition D( 3 )(tt„) is unlikely to 
hold for the considered GARCH(1,1) model. A more prominent decrease is 
observed within the proportions of anti-D( 4 )(u„) and anti-D^ 5 )(n n ). It will be 
seen in the simulation study that these proportions lead to a quite acceptable 
choice of values of k in the knew estimation procedure. 

Observe that from Proposition 12.51 we can also analyze D ( k \u n ) by evaluat¬ 
ing D( 2 )(w n ) within cycles {Z n } n >\. The respective plots are in Figures [3] and 
[I]and seem to corroborate the analysis above. 

The plots observation can give us some clue about D^ k \u n ) but does not 
allow to make a definite decision. We can always opt for higher values of k since, 
if D*- fc ) (u n ) holds then (u n ) holds for all s > k. However, a too large k for 
the cycles may diminish the precision of the new estimators, as will be pointed 
in the next section. 


4.1 Simulations 

In our study we consider 1000 replicates of simulated samples of size 1000 of 
each of the models referred previously: ARCauchy (p = —0.6), ARUnif (r = 2), 
MM (ao = 2/6, = 1/6,02 = 3/6), MAR (<f> = 0.5), Markov chain (a = 0.5) 

and GARCH(1,1) (A = 0.25, /3 = 0.7). We have calculated the values of the new 
estimator 9 given in (flHl) . as well as, the values of estimators 9 1 , 9 ML and 9 17 
based on the new indirect approach in ([171) . and estimator 9 ss based on (flT)l) . 
Although 9 ff and 9 FF are derived under a max-stable premise, we still apply 
them since, in practice, we are taking cycles {Z n } n >\ of maximums which, albeit 
crudely, can approach a max-stable behavior. We denote all these estimators 
as indirect. For comparison, we also consider the runs estimator (9 R ) and the 
intervals estimator (9 1 ) directly for 9x- In opposition to indirect estimators, we 
denote 9 R and 9 1 as direct. The root mean squared errors (rmse) and the abso¬ 
lute mean biases (abias) are given in Tables [7] and [3J for levels u n corresponding 
to the empirical quantiles 0.95, 0.975 and 0.99, respectively denoted, ( 70 . 95 , ( 70.975 
and ( 70 . 99 . For models MM, ARUnif, ARCauchy and MAR, which satisfy condi¬ 
tion D( 3 ^(m„), all the new estimators were based on the construction of cycles 
{Z n } n > 1 by taking k = 3. The direct runs estimator 9 R was also computed 
for run r = k (see Section [T|). These are reported in Table [7] The results for 
the Markov chain and GARCH (1,1) are given in Table [3l by considering that 
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Figure 1: From left to right and top to bottom, proportions of anti-D^ (it„) 
for ARCauchy, ARUnif, MM, MAR and Markov chain and anti-D^ 4 )(u„) of 
Markov chain, respectively, for r = 50 (full line) and t = 100 (dotted line), 
with k n = [(logn) 3 ]. 



Figure 2: From left to right and top to bottom, proportions of anti-D^ 3 ) (it ra ) 
with k n = [(logn) 3 ] and anti-D* 3 )(it ra ) to anti-D^ 5 )(it ra ) with k n = [(logn) 3 3 ] for 
GARCH(1,1), for r = 50 (full line) and r = 100 (dotted line). 
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Figure 3: From left to right and top to bottom, proportions of anti-D^ 2 )(u„) 
of cycles {Z n } n >i for ARCauchy, ARUnif, MM, MAR and Markov chain with 
k = 3 and Markov chain with k = 4, for r = 15 (full line) and r = 20 (dotted 
line), with k n = [(logn) 3 ]. 



Figure 4: From left to right and top to bottom, proportions of anti-D^ 2 )(u„) 
of cycles {Z n } n >i for GARCH(1,1) with k = 3 and k n = [(logn) 25 ] and 
GARCH(1,1) with k = 3,4,5 and k n = [(logn) 3 ], for r = 15 (full line) and 
r = 20 (dotted line). 
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D( fc )(u n ) holds with k = 4 in the first case and k = 5 in the second model (see 
Figures [T] to [3]). In what concerns the direct runs estimator 9 R , we choose a run 
r equal to k = 4 in the Markov chain model and run r equal to k = 5 in the 
GARCH(1,1) model. Ancona-Navarrete and Tawn ((T], 2000) considered r = 10 
for the runs estimator d R in the Markov chain model. Indeed, if we take r = 10 
in our simulations for this model, we obtain slightly lower rinse’s for this esti¬ 
mator. We have also considered r = 10 in the GARCH(1,1) model which led to 
an overall decreasing of 0.1 in the rmse’s for estimator d R . The presented choice 
of the values k for the indirect estimators leads to the best results among other 
values of k also tried in simulations but not reported. Indeed, if the models 
satisfy condition D( fc )(u„), the results by taking k + 1 are quite close but if we 
continue to increase k, they get worst. Observe that a too much large k means 
larger cycles {Z n } n >\ and thus some loss of information. 

The new approach presents good results, particularly with estimators 9 and 
9^ r . As expected, the upcrossings estimator is a competitor within our frame¬ 
work. The estimator 9 FF has also a good performance, except for the ARUnif 
model. In this case the results are better if we take k = 4, leading to a rmse 
ranging from 0.084 to 0.158. One reason is that the cycles {Z n } n > i with k = 4 
for this model may be more close of a max-stable behavior. We observe a sim¬ 
ilar situation with estimator 9 FF . It performs well except in model ARUnif 
where, for k = 4, we obtain a rmse of 0.077, as well as in model ARCauchy 
where k = 4 leads to a rmse of 0.063. The intervals estimator yields the largest 
errors and behaves better if applied indirectly in the case of the Markov chain 
and the GARCH(1,1). The indirect estimators 0 ML and 9 SS have a similar 
performance. The results tend to be better at lower thresholds. 

4.2 Application to financial data 

Log-returns of a financial time series usually present high volatility and cluster¬ 
ing of large values. Klar et al. ([15] , 2012) have analyzed DAX German stock 
market index time series and concluded that GARCH(1,1) is a good model to 
describe these data. In particular they considered the series of log-returns of 
DAX closing prices from 1991 to 1998 (see Figure[5]) and fitted a GARCH(1,1) 
model with autoregressive parameter A ~ 0.08, variance parameter /3 = 0.87 
and innovations 1 7 (after removing null log-returns). By the tabulated values of 
the extremal index of GARCH(1,1) models in Laurini and Tawn ([TB], 2012), the 
true value is around 0.3. In Table 0] we report the estimates, derived according 
to the conclusions of the simulations concerning the GARCH (1,1) model (see 
also the anti-D 1 ^ (m„) plots in Figures [ 6 ] and 0. Thus the direct runs estimator 
6 r was computed with run 5 and the indirect estimators (0, 0 17 , 9 1 , 0 ML , 9 ss , 
9 ff and 9 FF ) were calculated by considering cycles {Z n } n > 1 with k = 5. The 
closest values of 0.3 were obtained with quantile 0.95 in all cases, which is also 
in accordance with the simulation study. The indirect upcrossings estimator 9 U 
presents the nearest approximation, followed by the indirect intervals estimator 
9 1 and 9. We have also tried other values for k and found that, in this series, 
k = 6 leads to the best approximations of 0.3, with 9 = 0.34 and the remaining 
approximately 0.39, except for the intervals and the direct runs estimator where 
the estimates were 0.12 and 0.68. If we consider the direct runs estimator 9 R 
with run 10 (see Section 0) we obtain the estimate 0.48. 
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Table 2: The root mean squared error (rmse) and the absolute mean bias (abias) 
obtained for models MM, ARUnif, ARCauchy and MAR, by considering the 
empirical quantiles 0.95, 0.975 and 0.99, respectively, go. 95 , go .975 and go. 99 - 
The direct runs estimator 0 R is based on run r = 3. 



rmse 

abias 

MM 

go. 95 

go.975 

go.99 

go. 95 

go.975 

go.99 

qR 

0.055 

0.063 

0.095 

0.028 

0.002 

0.024 

e 1 

0.114 

0.152 

0.221 

0.061 

0.075 

0.137 

9 

0.057 

0.062 

0.095 

0.036 

0.007 

0.020 

& 

0.055 

0.077 

0.138 

0.013 

0.011 

0.058 

e 1 

0.141 

0.184 

0.268 

0.071 

0.105 

0.211 

qML 

0.063 

0.077 

0.138 

0.009 

0.008 

0.031 

qSS 

0.055 

0.084 

0.176 

0.000 

0.023 

0.112 

e FF 

0.032 

0.032 

0.055 

0.003 

0.004 

0.014 

9 ff * 

0.032 

0.003 

ARUnif 

go. 95 

go.975 

go.99 

go. 95 

go.975 

go.99 

qR 

0.063 

0.089 

0.138 

0.005 

0.011 

0.021 

e 1 

0.179 

0.130 

0.118 

0.200 

0.182 

0.202 

9 

0.003 

0.009 

0.019 

0.063 

0.089 

0.138 

9 F 

0.089 

0.118 

0.182 

0.011 

0.018 

0.039 

9 1 

0.130 

0.120 

0.145 

0.182 

0.195 

0.219 

qML 

0.015 

0.017 

0.022 

0.089 

0.122 

0.197 

9 SS 

0.020 

0.025 

0.088 

0.089 

0.118 

0.179 

9 ff 

0.335 

0.335 

0.335 

0.331 

0.331 

0.331 

9 ff * 

0.875 

0.861 

ARCauchy 

go. 95 

go.975 

go.99 

go. 95 

go.975 

go.99 

9 r 

0.077 

0.095 

0.152 

0.041 

0.013 

0.026 

9 1 

0.158 

0.182 

0.237 

0.095 

0.089 

0.132 

9 

0.084 

0.095 

0.152 

0.051 

0.019 

0.022 

& 

0.095 

0.130 

0.210 

0.018 

0.006 

0.068 

9 1 

0.179 

0.219 

0.286 

0.088 

0.112 

0.194 

qML 

0.095 

0.134 

0.219 

0.014 

0.006 

0.05 

qss 

0.089 

0.134 

0.219 

0.003 

0.026 

0.141 

9 ff 

0.084 

0.084 

0.084 

0.072 

0.074 

0.075 

9 ff * 

0.602 

0.595 

MAR 

go. 95 

go.975 

go.99 

go. 95 

go.975 

go.99 

9 r 

0.071 

0.095 

0.158 

0.005 

0.017 

0.058 

9 1 

0.134 

0.176 

0.261 

0.067 

0.087 

0.157 

9 

0.071 

0.094 

0.154 

0.026 

0.007 

0.051 

9° 

0.077 

0.114 

0.187 

0.009 

0.022 

0.072 

9 1 

0.145 

0.184 

0.251 

0.075 

0.105 

0.178 

qML 

0.077 

0.114 

0.192 

0.005 

0.021 

0.055 

qSS 

0.077 

0.114 

0.210 

0.002 

0.032 

0.126 

9 ff 

0.045 

0.055 

0.077 

0.003 

0.009 

0.023 

9 ff * 

0.032 

0.006 
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Table 3: The root mean squared error (rmse) and the absolute mean bias (abias) 
obtained for models Markov chain (MC) and GARCH(1,1), by considering the 
empirical quantiles 0.95, 0.975 and 0.99, respectively, < 70 . 95 ; 90.975 and 90 . 99 - The 
direct runs estimator d R is based on run r = 4 and r = 5 for, respectively, the 
MC and GARCH models. 



rmse 

abias 

MC 

90.95 

90.975 

90.99 

90.95 

90.975 

90.99 

qR 

0.084 

0.122 

0.202 

0.024 

0.053 

0.115 

e 1 

0.141 

0.184 

0.305 

0.082 

0.100 

0.202 

6 

0.071 

0.089 

0.141 

0.036 

0.013 

0.032 

e u 

0.071 

0.089 

0.141 

0.022 

0.005 

0.063 

e 1 

0.118 

0.148 

0.226 

0.032 

0.063 

0.136 

qML 

0.084 

0.110 

0.167 

0.021 

0.001 

0.063 

qSS 

0.077 

0.110 

0.187 

0.014 

0.017 

0.099 

e FF 

0.071 

0.077 

0.105 

0.053 

0.060 

0.078 

e FF * 

0.055 

0.050 

GARCH(1,1) 

90.95 

90.975 

90.99 

90.95 

90.975 

90.99 

e R 

0.148 

0.212 

0.295 

0.121 

0.175 

0.245 

e 1 

0.200 

0.221 

0.315 

0.130 

0.117 

0.215 

6 

0.110 

0.110 

0.141 

0.095 

0.075 

0.051 

9 ° 

0.105 

0.114 

0.152 

0.076 

0.057 

0.029 

e 1 

0.134 

0.126 

0.167 

0.085 

0.059 

0.003 

qML 

0.010 

0.110 

0.152 

0.073 

0.057 

0.027 

6 SS 

0.010 

0.110 

0.148 

0.071 

0.053 

0.002 

e FF 

0.063 

0.084 

0.134 

0.020 

0.019 

0.017 

e FF * 

0.045 

0.010 
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years 


Figure 5: Daily closing indexes (left) and daily log-returns (right) of DAX, from 
1991 to 1998, with 1786 observations (successive equal prices excluded). 


Table 4: Estimates of the extremal index of the DAX series at quantile 0.95. 
The direct runs estimator was derived with run 5. The indirect estimators (0, 
9 17 , 9 1 , 9 ML , 9 ss , 9 FF and 9 FF ) were obtained based on cycles {Z n } n >i with 
k = 5. 

9 r 9 1 9 9° 9 1 9 ml 9 ss 9 ff 9 ff * 

0.72 0.50 0.40 0.36 0.37 0.48 0.50 0.47 0.49 


5 Conclusions 

In this work we consider the estimation of the extremal index, an important 
dependence parameter within extreme values of stationary sequences. The new 
approach requires the validity of the local dependence condition D( fc )(u n ) of 
Chernick et al. ([3], 1991). The results are promising under a suitable choice 
for k and an empirical procedure was proposed for this evaluation. We also 
find that it is a useful tool for the well-known runs estimator, by guiding a first 
choice for the run. 
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Figure 6: From left to right, proportions of anti-D® (it ra ) to anti-D^ 6 ) (ii ra ) of 
the DAX series, with k n = [(logn) 2 ’ 5 ], for r = 15 (full line) and r = 20 (dotted 
line). 




d 


1 ' 
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